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Abstract 

In this paper we consider the production of top quark pairs close to threshold 
and compute the dependence on the Higgs boson mass to order aa s . This requires 
the evaluation of the matching coefficient of the vector current to two-loop level 
and the inclusion of Higgs mass dependent operators in the non-relativistic effective 
theory. For Higgs boson masses below 200 GeV moderate contributions to the top 
quark mass and the peak cross section are observed. We also provide additional 
information to the on-shell wave function renormalization which is relevant for the 
matching coefficient. 



PACS numbers: 14.65.Ha, 14.80.Bn 



1 Introduction 



One of the most important tasks of a future international linear collider (ILC) [1] is the 
precise measurement of the cross section for the production of top quark pairs close to 
threshold. The comparison with theoretical calculations will lead to a determination of the 
top quark mass with an unrivaled precision leading to an uncertainty below 100 MeV [2]. 
Furthermore, also the width of the top quark, the strong coupling and the Yukawa coupling 
between the top quark and Higgs boson can be extracted from such a measurement. The 
findings of this paper are important in this context. However, the final success of such an 
enterprise depends crucially on whether the theoretical precision can match the expected 
experimental one. The complete next-to-next-to-leading analysis of Ref. [3] has shown 
that it is important to obtain the third-order result within QCD. This constitutes a long- 
term calculation which has already been started (see, e.g., Refs. [4,5] for the most recent 
results). 

The naive scaling rule a ~ a 2 s) where a is the fine structure and a s is the strong cou- 
pling constant, shows that the next-to-next-to- next-to-leading terms are parametrically 
of the same order as the mixed corrections proportional to aa s . In this paper we take the 
first step to a complete order ota s result and compute the full dependence on the Higgs bo- 
son mass. On one hand this requires the evaluation of two-loop vertex corrections within 
the Standard Model (SM) involving a Higgs boson and a gluon. Furthermore, a new 
operator depending on the Higgs mass has to be considered on the effective theory side. 
As we will explicitly see, the two individual pieces are divergent, however, the physical 
quantities formed by the proper combination is finite. 

Effects from the Higgs boson mass have already been studied in Ref. [6] where a 
Yukawa-type potential has been used and the Schrodinger equation has been solved nu- 
merically in order to obtain the imaginary part of the Green function and finally the total 
cross section. In this paper we will consider the Higgs boson mass to be much larger 
than the soft scales involved in the process and include it in a systematic way in the 
non-relativistic effective theory. The importance of a systematic and careful treatment of 
electroweak corrections to the top quark pair production has been stressed in Ref. [7, 8] , 
where certain next-to-next-to-leading logarithmic electroweak effects associated to the 
instability of the top quark have been considered. 

The paper is organized as follows: In the next Section we consider the new Higgs 
mass dependent operator included into non-relativistic QCD (NRQCD) and evaluate the 
corrections to the energy level and the wave function at the origin. In Section |3] the 
Higgs mass dependent corrections to the matching coefficient are discussed and Section 0] 
contains a brief account of the phenomenological applications of our results. We conclude 
in Sectional In Appendix[X]we discuss the results for the matching coefficients induced by 
vertex corrections involving a Z boson and Appendix [B] contains additional information 
on the wave function renormalization constant which is relevant in connection to the 
matching coefficient. 
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2 Effective theory 



The standard theoretical framework for the threshold production of top quark pairs is 
given by NRQCD [9,10] which ensures the automatic resummation of terms like (a s /v) n 
in = 1,2,3,.. .) where v is the velocity of the produced quarks. NRQCD is constructed 
from QCD by integrating out the hard scales associated to the mass of the heavy quark, 
m. Thus it contains only degrees of freedom of the order mv and mv 2 . 

In this paper we consider next to the top quark, which takes over the role of the 
heavy quark, also the Higgs boson with mass Mh- Since Mh ^> m t v both mass scales are 
integrated out simultaneously. In Ref. [11] all operators which appear within QCD up to 
third order in perturbation theory have been classified. The only operator which has to 
be added to the Hamiltonian in Eq. (6) of Ref. [11] is given by 



anrm 2 



5Hh = ~s 2 M 2 M 2 ' (1) 

b W lvl W lvl H 

where sw is the sine of the Weinberg angle and My/ is the W boson mass. The cor- 
responding expression in coordinate space is proportional to the delta function. If we 
employ the counting rule a ~ a 2 it is easy to see that 6TCh gives contributions which are 
parametrically of the same order as the ones from the third-order QCD Hamiltonian. 

Let us in a first step evaluate the corrections to the energy levels induced by the 
operator 5H.H- Using first order perturbation theory we obtain 

«f - «ei*w»wf> - sf^gg^ e> 

with = —C F a 2 m t /(An 2 ) and is the Coulomb wave function, n is the principal 
quantum number. 

In order to compute the correction to the wave function at the origin the operator 
of Eq. (JTJ) has to be inserted into the standard formulae of non-relativistic perturbation 
theory 

<tyn(0) = - j d 3 f G c (0, f, E) 5H H V? if) , (3) 

where Go is the Coulomb Green function. Following Ref. [12] we split Gc into a con- 
tribution with zero, one and infinitely many gluon exchanges. Since only the one-gluon- 
exchange part is divergent it is convenient to perform the corresponding calculation in 
momentum space whereas the other contributions are evaluated in coordinate space. As 
a final result we obtain 



s W m W m h 



I f a s C F m t \ 3 
4 V H ) 8 



(4) 



with |^(0)| = C 3 F a 3 s m 3 t /(87m 3 ). The divergence in Eq. (J3J has been subtracted mini- 
mally which will also be done for the coefficient function considered in the next section. 

We want to mention that the formulae of this Section are adapted for the top quark. 
However, they also apply to other quark masses, in particular to the bottom system, by 
simply exchanging the top quark mass. 
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3 Higgs mass dependence of the matching coefficient 



The matching coefficient of the vector current j M = t'yH is defined by 

j k v = c^V x + o (J_\ , (5) 

where and x are two- component Pauli spinors for quark and anti-quark, respectively, 
and k = 1,2,3 denote the spacial components. Note that there is no contribution to our 
order from the time-component. 

For the practical computation of c v it is convenient to consider the tij vertex in the 
limit where for the photon energy, s, we have s ~ 4r??|. In this case we can establish the 
equation 

Z 2 T V = c v Z 2 Z~ 1 T v + ... , (6) 

where we have on the left- and right-hand side quantities of the full and effective theory, 
respectively. The latter are marked by a tilde and the ellipses denote terms suppressed 
by inverse powers of the top quark mass. r„ denotes the ti'y vertex corrections where 
it is understood that the couplings and masses are renormalized. The two-loop mixed 
correction to the on-shell wave function renormalization, Z 2 , has been computed recently 
in Ref. [13]. In particular, for the Higgs boson contributions both the exact result and 
the expansions in three kinematic regions have been provided. 

At this point it is convenient to apply the so-called threshold expansion [14, 15] to 
Eq. © which has the consequence that T v has to be evaluated for s = intf since all 
except the hard region cancel in Eq. (jBJ). Furthermore, on the right-hand side only tree 
contributions have to be considered. In particular we have Z 2 = 1. 

The one-loop corrections to c v are finite [16]. However, starting from two-loop order, 
the matching coefficient exhibits infra-red divergences which are compensated by ultra- 
violet divergences of the effective theory rendering physical quantities finite. In Eq. © 
the renormalization constant Z v which generates the anomalous dimension of j v takes 
over this part. Note that the vector current in the full theory does not get renormalized. 

It is convenient to introduce the perturbative decomposition of the matching coefficient 

Cv = i + ^c F cT^ l + (-) 2 c F cT^ 2 + ^c7 + ^c F c\ (7) 

where analog formulae also hold for Z 2 and T v . Since we consider in this paper only 
the contribution from the Higgs boson the corresponding quantities obtain an additional 
superscript H. 

The two-loop QCD corrections, c^ CD ' 2 , have been computed in Refs. [17, 18] (see also 
Ref. [19]). The complete SM corrections to one-loop order can be found in Ref. [16] where 
next to the vertex corrections also the box diagrams contributing to e + e~ — > ti have been 
considered. The latter, however, get no contributions from the Higgs boson. 

From the perturbative expansion of Eq. (JBJ) it is easy to see that the one-loop results 
c qcd,i an( j c ew are s j m piy given by the sum of the one-loop expressions for Z 2 and r„. 
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Whereas the individual pieces are still divergent the proper sum is finite. For convenience 
we repeat the QCD corrections and the Higgs mass dependent term of c^ w which are given 
by 



„QCD,1 



c 



H.ew 



-2 



m t 
M% 



3y 



^y 2 H 



2-9y% + I2y%_ ^ _ (-2 + 5^-6^) ^ f 



4m 



Hj 



where 



V4x 2 -1 



arctan \/Ax 2 — 1 for x > \ 



2x 2 111 



for x < h 



(9) 



2 ■ 



and = m t /M h . 

At two-loop order it is quite difficult to obtain a closed analytic result valid for all 
Higgs and top quark masses. However, it is possible to get compact formulae valid in 
various kinematical regions which - - when combined - - cover the whole Higgs mass 
range. This strategy has been successfully applied in Ref. [13] to the on-shell wave function 
renormalization (see also Appendix[BJ). Thus, let us consider c^ ,ew in the limits m t <C Mjj, 
m t pa Mh and m t 3> Mh where it is given by 



■'vAa 



Z vAb 



Ml 



y 2 H 



31 



+ y°H 



m t 

Ml 



144 24 
737 lllny| 



3 

16 



+ y 6 H 



307 5\ny 2 H 



480 



360 



6 



+ . . . 



71 I 1 

+ — — + VHAa 



6 8VS 



Tt 



+ Vhm 

1 



1 13tt 
36 108^ 



24 8x/3 



+ Vhm 



Vhm 



12 qVzJ 



+ 



57T 



288 ' 54^") 



+ 



rn 



M 2 W 



+ 



7T 



6 8^ 



VH,lb 



IT \ 



1 

24 + 8V3J 



Vh,u 



13 



19tt 



m t 
Ml 



+ 



7T 

4 yfl 
71 

720 ~ 



72 216^/ 
liay H 237T 1 



+ y%,ib 



J + Vh,u 
59 



8 24^/ 



1 

-7T + 



288 



23tt \ 
216y/3j 



2 

hi y H 
12 



96 y H 



+ 



7 3, 
1 — m y u 

48 8 



1 55tt 1 
S + 512^1 



1 

5£ 



(10) 



5 




Figure 1: c^ ,ew /(mf/M^) as a function of 1/ijh- The solid (black) line repre- 
sents the exact result. The approximations in the three regions are shown as 
dotted, dashed and dash-dotted lines where for the expansion around yn ~ 1 
two different approximations have been chosen. 



with y H ,\ a = (1 - 1/Vh) and VH,ib = (1 - Vh)- c„ f a w and c v £™ are two different rep- 
resentations of the same information which turn out to be useful for different values of 
Vh- 

In Fig. [U the exact result for c^ w /(mf/M^) (full black line, cf. Eq. (jHJ)) is shown 
together with the expansions from Eq. (JTUJ), c^ ew (dash-dotted), c^f a w (short dashed), 
c vib (dotted) and c^ w (long dashed). One can see that the approximations nicely 
cover the whole region of i/h- Note that the expansion around m t = Mn shows better 
convergence properties in case yu.ib is used as parameter. 

Let us draw the attention of the reader to the expansion terms for large top quark 
masses which starts with an enhancement factor vn? t / '(M^Mh) . Two out of the three 
powers in m t come from the Yukawa coupling between the Higgs boson and the top quark 
and the factor m t /Mh is the indication of the Coulomb singularity which would be present 
for a massless Higgs boson. The next-to-leading term is quadratic in mt accompanied by a 
logarithm in rnt/Mh- Both terms indicate that even for the leading terms it is not possible 
to nullify the Higgs boson mass and non-trivial integration regions have to be considered 
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Figure 2: Two-loop Feynman diagrams contributing to c. 



which makes the evaluation of the corresponding expansion terms quite tedious. This is 
particularly true for the two-loop order where due to the Coulomb divergence one will 
have a further factor rrit/Mh as compared to the one- loop term. Thus the expansion starts 
with a quartic top quark mass dependence. Furthermore, there are momentum regions 
which have \/m t / Mu as expansion parameter. 

On the other hand, as can be seen in Fig. ^ the result from this region is phenomeno- 
logically less important since for Higgs boson masses above 115 GeV there is perfect 
agreement between the exact result and the approximation one obtains for m t ~ and 
m t <C M h . Higgs boson masses below approximately 115 GeV are excluded by the direct 
search at the CERN Large Electron Positron Collider (LEP). For this reason we do not 
consider the limit mt ^> Mh at two-loop order. 

Let us now turn to two loops. The Feynman diagrams are shown in Fig. El They are 
generated with QGRAF [20]. The application of q2e and exp [21,22] identifies the topology 
of the individual diagrams and adopts the notation in order to match one of the following 
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functions 



J ± (ni,n2,n3,n4,n 5 ) = 

e 2e lE r d d kd d l 



(^/2)2 J (P + 2fcg)«i(/ 2 ±2/g)"2(A; 2 )n3((A;-/)2)n4(/2_ M 2)n 5 

H§{n 1 ,n 2 ,n 3 ,n4,n 5 ) = 

& 2e lE r d d kd d l 



( m d/2)2 J (k 2 + 2kq) n i{l 2 ±2lq) n *{k 2 ) n *{{k - I) 2 - M 2 ) n *{P) n s 
Y^(n 1 ,n 2 ,n 3 ,n 4z ,n 5 ) = 

e 2e7E f d d kd d l 



( m d/2)2 y (p + 2kq) n ^{l 2 ± 2lq) n *((k - I) 2 =F 2g(A; - l))^{k 2 )^{l 2 - M 2 )™* ' 
Z^(n 1 ,n 2 ,n 3 ,n i ,n 5 ) = 

e 2 ^ E f d d kd d l 
{m d l 2 ) 2 J (k 2 + 2kq) n ^(l 2 + 2lq)'^((k - I) 2 - 2q(k - Z))^(A; 2 - M 2 ) n4 (/ 2 )™ 5 ' ^ 

where d = 4 — 2e is the space-time dimension and the n« are integer indices. The choice 
of the five-line integrals of Eq. (jllj) is possible due to the special kinematics we have at 
hand. Note that J + , H + and Y~ correspond to self energies whereas J~ , H~ , Y + and 
Z + to vertex diagrams. 

In a next step we use the program AIR [23] in order to reduce the integrals to about 30 
master integrals. They range from very simple two-point expressions up to complicated 
two-scale integrals with five lines. At this point an asymptotic expansion in the various 
kinematical regions is performed. From the one-loop result (cf. Fig. ^) and from the 
considerations in the context of the on-shell wave function renormalization (see Ref. [13] 
and Appendix [Bj) one can see that for the phenomenological interesting Higgs boson 
masses it is sufficient to have the expansion around m t « Mh and the one for m t <C 
Mh at hand. Thus the master integrals are expanded in these two limits. In the case 
m t ~ Mh this reduces to a simple Taylor expansion whereas for m t <C Mh the well- 
established hard-mass procedure [15] is applied. The latter is actually automated in the 
program exp [21,22]. Hence as an independent check we applied with the help of exp the 
hard-mass procedure to each diagram which immediately leads to simpler expressions and 
avoids the reduction to the master integrals. In this way the result for m t <C Mh could be 
checked. We want to mention that the calculation has been performed for general QCD 
gauge parameter, £. The cancellation of £ in the final result serves as a welcome check 
for our calculation. 

As already mentioned in the Introduction, even after the proper combination of Z 2 
and is formed, as prescribed by Eq. (|SjL there remains an infra-red divergence which 
is plugged into Z v . We subtract this divergence in the MS scheme and introduce the 
anomalous dimension j v = which is given by 

_ OLOi s tt 2 m\ 

lv - ~^ Cf ^mimi- (12) 
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In the regions m t <C Mh and m t ~ Mh three, respectively, six expansion terms have 
been evaluated. Our results read 
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(13) 



where m t is the on-shell mass and (3 ~ 1.20206 is Riemann's zeta-function. The ln(m 2 //x 2 ) 
term reminds on the divergence which has been subtracted minimally. Since in the limit 
m t <C Mh some coefficients of the e-expansion of the master integrals could only be 
computed numerically the results for the finite parts of c^™ lx and c^™ lx are presented 
in numerical form. We also want to stress again that c^™ lx and c v '^ x contain the same 
information. However, expressed in terms of yn,\b the convergence properties are much 
better. 

In Fig. El the results are shown where the same notation as for the one- loop result has 
been adopted and \i = m t has been chosen. Next to the highest expansion terms we also 
include the lower-order terms as thin lines which nicely demonstrates the convergence 
properties in the individual regions. Comparing the two parameterizations of the expan- 
sion around M t = Mh one observes that the one in terms of yH,u demonstrates a much 
better convergence behaviour: the thin dots in Fig. |3] are practically indistinguishable 
from the fat ones. This can also be seen in Eq. (J13)) where the coefficients become quite 
small starting from the third term which is not the case for c^™ lx . Thus, the combination 



of c^™ lx and c^Q mix provide a very good approximation to for Higgs boson masses 
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Figure 3: c™ lx as a function of 1/vh- The approximations are shown as dotted, 
dashed and dash-dotted lines where for the expansion around yu ~ 1 two 
different parameterizations have been chosen. Lower-order terms are plotted 
using thin lines. For the renormalization scale \i = m t has been adopted. 



4 Phenomenological application 

Let us in a first step discuss the effect on the top quark mass. The connection between 
the position of the peak in the threshold cross section, E res , and the top quark mass is 
given by 

E ves = 2m t + Ef- t - + 5 Tt E Tes , (14) 

where Ef' ' is the perturbative part of the ground state energy and 5 Ft E ms = 100 ± 
10 MeV [24] takes into account the effect of the finite width, the higher order resonances 
and the continuum. Non-perturbative effects are negligible for the top quark system. Ef' ' 
up to third order in QCD has been computed in Ref. [24]. The contribution from the 
Higgs boson is given in Eq. (J2J) where ot s has to be evaluated at the typical soft scale given 
by jji s = Ci?a s (/^ s )mt ~ 30 GeV. 

According to Eq. ()14J) a shift in E\' ' can directly be translated into a shift in the 
top quark mass, 5m t . Choosing for the input values m t = 175 GeV, M w = 80.41 GeV, 
s 2 w = 0.23, a = 1/127 [25], a s (M z ) = 0.118, M H = 120 GeV and /i s = 30 GeV we observe 
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for 8m t « -E? /2 « 26 MeV. This reduces to 9(1) MeV for M H = 200(500) GeV. If 
M# = 120 GeV is adopted and fi s = 15(60) GeV is chosen one obtain 8m t rs 38(18) MeV. 
Thus, for light Higgs masses relatively large effects are observed whereas for larger M# 
the numerical effect on rrit is small. Note that our findings are in agreement with the 
more qualitative analysis of Ref . [6] . 

As a second application it is interesting to consider the effect of the new corrections 
to the normalized cross section R = a(e + e~ — > tF)/a(e + e~ — > /U + /i~) at the resonance 
energy which is dominated by the contribution from the would-be toponium ground-state. 
It is of the form 



p _ pLO 2/ \ l^i(Q)! 2 | n ^ 



with R\° = 67[N c Qf\ipf(0)\ 2 / (mfT t ). The ellipses in Eq. (|TH|) denote contributions from 
higher order operators. Note that the divergences in c v and ipi{0) cancel in the combi- 
nation of Eq. (|15p. The most advanced evaluation of R\ is provided in Ref. [4] where 
all logarithmically enhanced third-order corrections and the ones proportional to 0q are 
included. Let us for completeness repeat the final numerical result which is given by 

R 1 « ^ o (l- 0.244 NLO + 0.449 NN LO- 0.277 N;iL o' + 5g ) +4 2) ...), (16) 

where /x s = 30 GeV has been chosen. The prime reminds that the N 3 LO corrections 
are not complete and 5$ and 8jf parameterize the one- and two-loop corrections due to 
the Higgs boson considered in the present paper. Using Mh = 120/200/500 GeV and 
Us = 30 GeV we obtain <?g } = 0.067/0.034/0.009 and = 0.036/0.011/0.0002. Thus, 
moderate effects are observed for Higgs boson masses below about 200 GeV. However, in 
contrast to the pure QCD effects of Eq. (fTUj) the convergence seems to be much better as 
8 $ it is substantially smaller than 8ff. 

In principle it is possible to apply the formulae of this Section also to the bottom 
system. However, due to the suppression factor ml/Mfj the numerical effect is very small 
for Higgs boson masses above 100 GeV. 



5 Conclusions 

In this paper we take an important step towards the evaluation of the mixed 
QCD/electroweak corrections to the threshold production of top quark pairs, namely, 
the complete Higgs boson mass dependence is computed. In particular, the two-loop 
corrections to the matching coefficients are evaluated and the corresponding operator is 
introduced in the effective theory. Moderate numerical effects on the top quark mass and 
the peak cross section are observed for Higgs boson masses below 200 GeV. Considering 
the anticipated precision of an ILC it is certainly necessary to take into account these 
corrections. 
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A Matching coefficients for the Z boson exchange di- 
agrams 

In the following we present the results for the matching coefficients induced by vertex 
corrections involving the Z boson and the corresponding Goldstone boson, respectively. 
The calculation is similar to the diagrams involving a Higgs boson which is presented in 
the main part of the paper. Note that there are also box diagrams involving the Z boson 
which are not considered here. The one-loop result is given by [16] 
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with y z = m t /M z , a t = l/(2s w c w ), v t = (1/2 - 4s^ v /3)/(s w c w ) and c w = y/T- s\ 
Mw/Mz- The function ^(x) can be found in Eq. @. It is straightforward to obtain 
the expansion in the various kinematical regions. Thus we refrain from listing them 
explicitly. In Fig. Ufa) the exact result is plotted together with the expansions where in 
each region the same depth is chosen as in Eq. (fTU|) . The vertical dashed line marks the 
phenomenological result for yz which is nicely approximated both from the m t 3> Mz 
and the m t w Mz region. However, it should be noted that both parameterizations for 
the m t = Mh expansion become instable below l/yz ~ 0.6. 

At two-loop order we consider the diagrams of Fig.[2]where the Higgs boson is replaced 
by the Z boson. The reduction of the occuring integrals is in close analogy and in 
fact the same set of master integrals is necessary. Although the region m t <ti Mz is 
phenomenologically not relevant we nevertheless present the results since it constitutes 
an important cross checks for the other kinematical limits. Furthermore, it is possible 
to obtain from this limit the result for the case of the bottom quark. We adopt the 
notation from Eq. ((7j) and denote the contribution from the Z boson to c™ lx by cf' mix . 
The additional subscript "0", "la" or "lb" reminds on the kinematical region considered. 
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Figure 4: (a) cf ,ew and (b) cf ' mix as a function of 1/yz- The same notation as 
in Figs. H and El is adopted. The vertical line indicates the physical value of 

l/yz = M z /m t . 



Our results read 



c 



^Z,mix 



7r 2 , m? 



/ 2 i 2\ 2 " l ""t 

(u t +a t ) s^yln — 



, 2 2 



587 131tt 2 
" 576 



vr 2 ln2 21 



28 



11 



C 3 + ^-\ny z -^-\n 2 y 2 z )y z 



432 



16 



27 



72 



+ 



+ 



403 205vr 2 



+ 



2299 2 5 2 2 i i 
+ W ln ^-18 ln 



64 576 

21091951 3951tt 2 28301 



, 2 2 

-\-a t s w 



864000 
281 



+ 



432 



+ 



2560 
85tt 2 



+ 



469 



In yi + — In 2 yl I y 7 
7200 z 160 z ' z 



7r 2 ln2 21 



29 



11 



-Z( 3 + --lny z -^-ln 2 y 2 z )y z 



+ 



+ 



709 llvr 2 
+ + 



192 
521 



79 



16 



ln 2/z- i^^vl ) Vz 



54 



72 



576 432 288 72 
/260687 16013vr 2 32453 



32000 



+ 



69120 21600 



5479 
1440 



\^y 2 z)y z + ... 



13 



^Z,mix 



2 9 

7r . mi 



2/Z,la 



, 2 2 



- 4.564 - A.699y ZM - 5.009y 2 ZM - 5.277y% la - 5.502y% la 



-5.69%* + 



.22 
' a t S W 



- 1.324- 1.504^- 1.7402/1^ 



-1.930y| la - 2.083y| ilo - 2.212y| tla + 



Z,mix 



7T^ 771^ 

K 2 + a?) s^y In -|- (1 - j/ z>16 ) + v 2 t s 



2„2 



- 4.564 + 4.699y Zi i fe 



-0.311y| fl6 - 0.043y| lb + 0.001^ 16 + 0.011y| fl6 + 



, 2 2 



- 1.324 + 1.504^, 16 - 0.236y| il6 - 0.047y| )lft - O.Ollyf^ 



-0.001y|, 16 + 



where the divergence has been subtracted in a minimal way and yz = m t /Mz, yz,ia = 
(1 — 1/yz) and yz,u = (1 — y%)- I* 1 Fig- HP 3 ) the results of Eq. (fTHj) are shown as a function 
in 1/ . Similarly to the two-loop case one observes a rapid convergence for c v '^ 1X whereas 
the magnitude of the coefficients of the higher order terms in cf'^ lx even increase. This 
gives us quite some confidence that c z v '^ lx evaluated at the physical scale provides a very 
good approximation to the unknown exact result. Using 1/yz = 91.19/175 ~ 0.52 we get 



„Z,mix 



-3.1 ±0.3 



(19) 



where we assign a generous uncertainty of 10% which results from the comparison of 
lower-order approximations. The accuracy of the result for c^' mix as given in Eq. (J 19)) is 
more than sufficient as far as the forseen precision of the measurement of the threshold 
top quark production is concerned. 



B Addendum to Ref. [13] 

In this Appendix we want to provide an update of the results of Ref. [13]. In Fig. 4 of that 
reference the Higgs boson mass dependence of the on-shell wave function renormalization 
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(c) 

Figure 5: (a) One-loop, (b) 1/e pole and (c) constant part of the two-loop 
corrections to Z^'°^ as a function of 1/yH — Mjf/mt. The solid (black) line 
represents the exact result. The approximations in the three regions are shown 
as dotted, dashed and dash-dotted lines where for the expansion around ys « 1 
two different parameterizations (short dashes and dots) have been chosen. 



constant to one- and two-loop order has been plotted as a function of 1/uh — Mjj/mt- 
Next to the exact results also the ones for m t <C Mh, tth ~ Mh and m t 3> Mh have been 
shown and good agreement over almost the whole range in yn has been found. However, 
there was a small gap around 1/yn ~ 2 which was not covered very well. This deficit 
is removed in Fig. El where the result obtained in the limit m t m Mh is plotted both in 
terms of yn,i a = (1 ~ ^-/Vh) (dashed) and yn,n — (1 ~ Uh) (dotted). One can see that 
the former is valid down to quite small values of 1/yn whereas the validity of the latter 
expansion reaches up to larger values of 1/yn leading to a significant overlap with the 
expansions for large and small values of Thus, in the whole yu range the simple 
expansions approximate the exact result to a very high precision. 

Let us also mention that there is a misprint in the definition of v% in Ref. [13] after 
Eq. (17): it is too small by a factor two. Furthermore, the analytic expression for Y e in 
Eq. (22) has to be multiplied by (—1) and in Eq. (27) a factor is missing in the 



15 



terms containing a s /n. 



References 



[i] 

[2] 
[3] 
[4] 

[5 

[6] 
[7] 



[9 
[10 

[11 
[12 

[13 
[14 
[15 

[16 

[17; 



see, e.g., http://linearcollider.org/ 

M. Martinez and R. Miquel, Eur. Phys. J. C 27 (2003) 49. 

A. H. Hoang et al, Eur. Phys. J. directC 2 (2000) 1 ar Xiv:hep-ph/0001286| . 

A. A. Penin, V. A. Sm irnov and M. Steinhauser, Nucl. Phys. B 716 (2005) 303 
|arXiv:hep-ph/0501042| . 

M. Beneke, Y. Kiy o and K. Schuller, Nucl. Phys. B 714 (2005) 67 
jarXiv:hep-ph/0501289| . 

M. J. Strassler and M. E. Peskin, Phys. Rev. D 43 (1991) 1500. 

A. H. Hoang and C. J. Reisser, Phys. Rev. D 71 (2005) 074022 
|arXiv:hep-ph/0412258| . 

A. H. Hoang and C. J. Reisser, arXiv:hep-ph/0604104. 

W. E. Caswell and G. P. Lepage, Phys. Lett. B 167 (1986) 437. 

G. T. Bodwin, E. Braaten and G. P. Lepage, Phys. Rev. D 51 (1995) 1125 [Erratum- 
ibid. D 55 (1997) 5853] |arXiv:hep-ph/9407339| . 

B. A. Kniehl, A. A. Penin, V. A. Smirnov and M. Steinhauser, Nucl. Phys. B 635 

(2002) 357 |arXiv:hep-ph/0203166| . 

B. A. Kniehl, A. A. Penin, M. Steinhauser and V. A. Smirnov, Phys. Rev. Lett. 90 

(2003) 212001 |arXiv:hep-ph/0210161| . 

D. Eiras and M. Steinhauser, JHEP 0602 (2006) 010 | |arXiv:hep-ph/0512099| . 



M. Beneke and V. A. Smirnov, Nucl. Phys. B 522 (1998) 321 |arXiv:he p-ph/9 711391| . 

V.A. Smirnov, Applied Asymptotic Expansions in Momenta and Masses, Springer- 
Verlag, Berlin- Heidelberg, 2001. 

R. J. Guth and J. H. Kiihn, Nucl. Phys. B 368 (1992) 38. 

A. Czarnecki and K. Melnikov, Phys. Rev. Lett. 80 (1998) 2531 
|arXiv:hep-ph/9712222| . 



[18] M. Beneke, A. Si gner and V. A. Smirnov, Phys. Rev. Lett. 80 (1998) 2535 
|arXiv:hep-ph/97123"02| . 



16 



[19] B. A. Kniehl, A. Onishchenko, J. H. Piclum and M. Steinhauser, 
|arXiv:hep-ph/0604072| Phys. Lett. B (in print). 

[20] P. Nogueira, J. Comput. Phys. 105 (1993) 279. 

[21] R. Harlander, T. Seidensticker and M. Steinhauser, Phys. Lett. B 426 (1998) 125 
|hep-ph/9712228| . 

[22] T. Seidensticker, |hep-ph/9905298| 

[23] C. Anastasiou and A. Lazopoulos, JHEP 0407 (2004) 046 | arXi v:hep-ph/0404258| . 

[24] A. A. Penin and M. Steinhauser, Phys. Lett. B 538 (2002) 335 
|arXiv:hep-ph/0204290| . 

[25] F. Jegerlehner, |arXiv:hep-ph/0105283| 



17 



